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Abstract 

The first measurements of the polarized structure function (t^x' for the reaction p{e,e' K~^)A. 
in the nucleon resonance region are reported. Measurements are included from threshold up to 
VF=2.05 GeV for central values of Q'^ of 0.65 and 1.00 GeV^, and nearly the entire kaon center- 
of-mass angular range. (Jlt' is the imaginary part of the longitudinal-transverse response and is 
expected to be sensitive to interferences between competing intermediate s-channel resonances, as 
well as resonant and non-resonant processes. The results for a^x' are comparable in magnitude to 
previously reported results from CLAS for glt^ the real part of the same response. An intriguing 
sign change in a^T' is observed in the high data at ~ 1.9 GeV. Comparisons to several 
existing model predictions are shown. 

PACS numbers: 13.40.-f, 13.60.Rj, 13.88. +e, 14.20.Jn, 14.40.Aq 
Keywords: kaon electroproduction, polarization, structure functions 
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I. INTRODUCTION 



The study of the electromagnetic production of strange quarks in the resonance region 
plays an important role in understanding the strong interaction. The p{e, e'K~^)A reaction 
involves the production of the strange particles A{uds) and K~^{us) in the final state via 
strange quark-pair (ss) creation. The fundamental theory for the description of the dy- 
namics of quarks and gluons is known as quantum chromo dynamics (QCD). However, while 
numerical approaches to QCD in the medium-energy regime do exist, neither perturbative 
QCD nor lattice QCD can presently predict hadron properties seen in this type of reaction. 
In the non-perturbative regime of nucleon resonance physics, the consequence is that the 
interpretation of dynamical hadronic processes still hinges to a significant degree on mod- 
els containing some phenomenological ingredients. Various quark models (see for example 
Refs. IJ]) predict a large number of non-strange baryons that can decay into a 

strange baryon and a strange meson, as well as Ntt/Ntttt final states. While many of these 
excited states have been observed in pion production data, a large number are "missing". 
The higher threshold for K~^A final states kinematically favors production of the missing 
resonances with masses near 2 GeV. Studies of different final states, such as the associ- 
ated production of strangeness, can provide complementary information on the contributing 
amplitudes. 

In the absence of direct QCD predictions, effective models must be employed. Utilizing 
these models by means of fitting them to the available experimental data - cross sections 
and polarization observables - or comparing the data to the model predictions, can provide 
information on the reaction dynamics. In addition, these comparisons can provide impor- 
tant qualitative and quantitative information on the contributing resonant and non-resonant 
terms in the s, t, and u reaction channels (see Fig. [1]). The development of these theoretical 
models has been highly based on the availability of the experimental data. Precise measure- 
ments of cross sections and polarization observables are crucial for the refinement of these 
models and for the search for missing resonances. 

In this paper, we report the first-ever measurements of the longitudinal-transverse polar- 
ized structure function, aiT', for the p{e, e'K~^)A reaction in the resonance region, using the 
CEBAF Large Acceptance Spectrometer (CLAS) in Hall B of Jefferson Lab. This observable 
provides complementary information to the cr^y structure function reported in Ref. j^, as 
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s-channel t-channel u-channel 

FIG. 1: Feynman diagrams representing s-channel nucleon {p, N* , A*) exchange (left), t-channel 
kaon {K,K*) exchange (middle), and u-channel hyperon {Y,Y*) exchange (right) that contribute 
to the reaction models. The vertex labels qmbb represent the strong coupling constants. 



will be discussed. Thus, these new data provide another constraint on model parameters, 
and therefore, provide additional important information in understanding the process of 
electromagnetic production of strangeness. 

There is a growing body of high-quality data on the electromagnetic production of strange 
hadrons. Recently published data using electron beams exist on the separation of the lon- 
gitudinal and transverse structure functions, ai and ar, from Hall C of Jefferson Lab for 



both K~^A and K~^T,^ final states at 
of-mass scattering angle of 6'^ = 0° 



/V=1M GeV, for up to 2.0 GeV^, at a kaon center- 
6|, lZI] . The CLAS Collaboration has recently produced 
results in which unpolarized cross sections and interference structure functions {(Ttt and 
air) have been measured for K~^A and K~^J]^ final states over a wide kinematic range with 
up to 2.6 GeV^, W up to 2.4 GeV, and nearly complete angular coverage in the center-of- 
mass frame These results include the first-ever separation of and o"t at angles other 
than 6*1^' = 0°. The same set of data has been analyzed to extract the polarization transfer 
from the virtual photon to the produced A hyperon 8|] and to extract the ratio of o"l/o"t 
at 6*1^ = 0° for the K~^A final state 9|. Older electroproduction data from various labs also 



exist 



[lO 
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but with much larger uncertainties and much smaller kinematic 



coverage. 

Complementary data from photoproduction are also available. The SAPHIR collabora- 
tion has published total and differential cross section data for photoproduction of K~^A and 



K'^T,^ final states with photon energies up to 2 GeV 15,U6|. CLAS has provided extensive 

nn r n 

differential cross sections [17|, Il8[ , along with recoil [17[ and transferred polarization [19|] data 
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for the same final states in similar kinematics. Finally, the LEPS collaboration has measured 
differential cross sections and polarized beam asymmetries with a linearly polarized photon 
beam for energies up to 2.4 GeV at forward angles 2^, 21]. 

The SAPHIR cross section data show an interesting resonance-like structure in the h. 
final state around W=\.9 GeV. A similar structure has been seen in the unpolarized elec- 
troproduction cross section data j^, as well as in the photoproduction measurements of 
CLAS Q,Q. Within the isobar model of Mart and Bennhold |22l . |23| ] , that structure was 
interpreted as a Di3(1895) resonance, which had been predicted by several quark models 
(e.g. Ref. ^), but not well established. However, the isobar model of Saghai 24 1 found that 
the cross section data could be satisfactorily described without the need for including any 
new s-channel resonances by including higher-spin u-channel exchange terms. The need to 
include the missing Z)i3(1895) state, however, was supported by the new Regge plus reso- 



nance model from Ghent 



25l | that compared their model to a broad set of cross section and 



polarization observables from the available photo- and electroproduction data. 

The organization of this paper includes an overview of the relevant formalism in SectionllTl 
a description of the theoretical models used to compare against the data in Section UTTl details 
on the experiment and data analysis in Section HVt and a presentation and discussion of the 
results in Section |Vl The conclusions are given in Section IVll 



II. FORMALISM 

A schematic diagram of K^K electroproduction off a fixed hydrogen target is shown in 
Fig. [21 The angle between the incident and scattered electron is 6'e, while the angle between 
the electron scattering plane and hadron production plane is defined as </>. In the one-photon 
exchange approximation, the interaction between the incident electron beam and the target 
proton is mediated by a virtual photon, 7*. The virtual photon four momentum is obtained 
from the difference between the four momenta of the incident, e = (-E,Pe), and scattered 
electrons, e' = (i?',Pg), as: 

g = e - e' = (z/,q). (1) 
The four momentum transfer squared, Q^, is an invariant quantity defined as: 

Q2 = -q^ = -{u^- q2) = AEE' sin2(e,,/2), (2) 
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FIG. 2: (Color online) Kinematic diagram for kaon-hyperon (KY) electroproduction. 

where u = E — E' is the energy transfer and 9^1 is the electron scattering angle in the lab 
frame. The invariant mass W of the intermediate hadronic state is defined as: 

W'^ = s = M'^ + 2Mv - (3) 



where M is the mass of the proton 



target. 



Following the notation of Refs. 26|, |27l], the differential cross section for KY electropro- 



duction in the center-of-mass frame is given by: 



dE'dn'jn*j^ dQj^ 
where F is the virtual photon flux given by: 

a E' W^-M^ 



da ^ 



(5) 



47r2 EM Q2 

Here is the virtual photon cross section and e is the virtual photon transverse polariza- 
tion component defined as: 

The cross section for the electromagnetic interaction of a relativistic electron beam with 
a hadron target is obtained by calculating the transition probability of the process 28| . The 
cross section can be written in the form of a contraction between leptonic and hadronic 
tensors that contain the electron and hadron variables separately. In general, the lepton 
tensor can be written in terms of a density matrix of virtual photon polarization that con- 
tains a symmetric helicity- independent part and an anti-symmetric helicity-dependent part. 
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The anti-symmetric part contributes to the cross section only when the hadron tensor also 
contains an anti-symmetric part. This is the case when scattering a polarized electron off 
of an unpolarized target with the detection of the final state hadron in coincidence with the 
scattered electron. The anti-symmetric part vanishes for the case of unpolarized electrons. 

For a polarized electron beam with helicity h and no target or recoil polarizations, the 
virtual photon cross section can be written as: 

da 



dn*j^ 



ax + eaL + eaxr cos 20 + a/ e(l + e)aLT cos + /i ye^T— e)(TLT' sin 0, (7) 



where ai are the structure functions that measure the response of the hadronic system 
and i = T, L, LT, TT, and LT' represent the transverse, longitudinal, and interference 
structure functions. The structure functions are, in general, functions of Q^, W, and 9}^ 
only. Note that the convention employed here for the differential cross section is not used 
by all authors 29 1. 



For the case of an unpolarized electron beam, Eq.Q reduces to the unpolarized cross 
section, ctq: 

da 



dn*K 



ao = aT + (^o-L + eo"rT cos 20 + iy/e(l + €)aLT cos 0. (8) 



The electron polarization therefore produces a fifth structure function that is related to the 
beam helicity asymmetry via: 



dcr ' da / — r r- 



Alt' = 7+ 7- = ^ (9) 

dn}^ ~^ dn*^. 

The ± superscripts on correspond to the electron helicity states of h = ±1. Clearly, a it' 
can only be observed when the outgoing hadron is detected out of the electron scattering 
plane (0 ^ 0) and can be separated by flipping the electron helicity. 

The structure functions are defined in terms of the independent elements of the hadron 



tensor in the center-of-mass frame, W^y 28l |: 



(Tl oc 1^00' 

ar oc (W^n + 

(Ttt oc VTi'.i, (10) 

(Tlt oc 

aLT' cx 



where the indices A, A' = for the longitudinal component and A, A' = ±1 for the two 
transverse components. In contrast to the case of real photons, where there is only the 
purely transverse response, virtual photons allow longitudinal, transverse-transverse, and 
longitudinal-transverse interference terms to occur. 

The polarized structure function ctlt' is intrinsically different from the four structure 
functions of the unpolarized cross section. As seen by Eqs. fllOl) . this term is generated by 
the imaginary part of terms involving the interference between longitudinal and transverse 
components of the hadronic and leptonic currents. This is in contrast to (Jlti which is 
generated by the real part of the same interference. aiT' is non-vanishing only if the hadronic 
tensor is anti-symmetric, which will occur in the presence of final state interaction (FSI) 
(or rescattering) effects, interferences between multiple resonances, interferences between 
resonant and non-resonant processes, or even between non-resonant process alone. On the 
other hand, ar^T' could be non-zero even when aLT (which is not expected to be sensitive 



to FSI effects 



281]) is zero. It provides a means of measuring the contributions of small 
resonance channels that are often too weak to be observed directly in the unpolarized cross 
sections. Furthermore, when the reaction proceeds through a channel in which a single 
amplitude dominates, the longitudinal-transverse response will be real and aj^T' vanishes. 
Both (Jlt and ctlt' are necessary to fully unravel the longitudinal-transverse response of the 
K^h. electroproduction reaction. 

III. THEORETICAL MODELS 

With the recently available data from the photo- and electroproduction of KY final 
states from CLAS and elsewhere, there have been renewed efforts on the development of 
theoretical models. The majority of these are single-channel models that represent tree-level 
calculations, where the amplitude is constructed from the lowest-order Feynman diagrams 
(see Ref. 25| and references therein). More recent work has moved beyond the single- channel 



approach with the development of coupled- channels models 3d, l3l|, 122, l33| or by fitting 




simultaneously to multiple, independent reaction channels |34l. I35l|. However, as a combined 
coupled-channels analysis of the photo- and electroproduction reactions is not yet available, a 
tree-level approach currently represents the best possibility of studying both reactions within 
the same framework. While most of the recent theoretical analyses have focused solely on 
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the available photoproduction data, it has been shown that electroproduction observables 



25|. 



can yield important complementary insights to improve and constrain theory 

At the medium energies used in this experiment, perturbative QCD is not capable of 
providing any predictions for the differential cross sections or structure functions for kaon 
electroproduction. In this work, the results are compared against three different model 
approaches. The first is a traditional hadrodynamic (resonance) model, the second is based 
on a Reggeon-exchange model, and the third is a hybrid Regge plus resonance approach. 

In the hadrodynamic model approach, the strong interaction is modeled by an effective 
Lagrangian, which is constructed from tree-level Born and extended Born terms for inter- 
mediate states exchanged in the s, t, and u reaction channels (see Fig. [T]). Each resonance 
has its own strong coupling constants and strong decay widths. A complete description of 
the physics processes requires taking into account all possible channels that could couple to 
the initial and final state measured, but the advantages of the tree-level approach include 
the ability to limit complexity and to identify the dominant trends. In the one-channel, 
tree-level approach, several dozen parameters must be fixed by fitting to the data, since 
they are poorly known and not constrained from other sources. 

The hadrodynamic model employed in this work was developed by Mart and Bennhold 



23|, |36|] (referred to here as MB). In this model, the coupling strengths have been determined 



mainly by fits to existing 7p — * K~^Y data (with some older electroproduction data included), 
leaving the coupling constants as free parameters (constrained loosely by SU(3) symmetry 
requirements) . It employs phenomenological form factors to account for the extension of the 
point-like interactions at the hadronic vertices. This model has been compared against the 



existing photoproduction data from SAPHIR 15|, ll6| and CLAS 



17 



and provides a fair 



description of those results. The model parameters are not based on fits to any CLAS data. 
The specific resonances included in this model are the S'ii(1650), Pii(1710), Pi3(1720), and 
^13(1895) A^* states in the s-channel, and the i^'*(892) and i^j'(1270) in the t-channel. 

The data are also compared to the Reggeon-exchange model from Guidal, Laget, and 
Vanderhaeghen jsi] (referred to here as GLV). This calculation includes no baryon resonance 
terms at all. Instead, it is based only on gauge-invariant t-channel K and K* Regge- 
trajectory exchange. It therefore provides a complementary basis for studying the underlying 
dynamics of strangeness production. It is important to note that the Regge approach has 
far fewer parameters compared to the hadrodynamic models. These include the K and K* 
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form factors (assumed to be of a monopole form) and the coupling constants Qkyn and 
Qk'YN (taken from photoproduction studies). 

The GLV model was fit to higher-energy photoproduction data where there is little doubt 
of the dominance of these kaon exchanges, and extrapolated down to JLab energies. An 
important feature of this model is the way gauge invariance is achieved for the K and 
K* t-channel exchanges by Reggeizing the s-channel nucleon pole contribution in the same 



manner as the t-channel diagrams 



371 ]. Due to gauge invariance, the t-channel exchanges 



and s-channel nucleon pole terms are inseparable and are treated on the same footing. 
They are Reggeized in the same way and multiplied by the same electromagnetic form 
factor. No counter terms need to be introduced to restore gauge invariance as is done in the 
hadro dynamic approach. 

The final model included in this work was developed by the University of Ghent group 25| , 



and is based on a tree-level effective field model for A and TP photoproduction from the pro- 
ton. It differs from traditional isobar approaches in its description of the non-resonant 
diagrams, which involve the exchange of K and K* Regge trajectories. A selection of 
s-channel resonances are then added to this background. This "Regge plus resonance" 
(referred to here as RPR) approach has the advantage that the background diagrams con- 
tain only a few parameters that are constrained by high-energy data where the t-channel 
processes dominate. Furthermore, the use of Regge propagators eliminates the need to intro- 
duce strong form factors in the background terms, thus avoiding the gauge-invariance issues 
associated with the traditional effective Lagrangian models. In addition to the kaonic trajec- 
tories to model the t-channel background, the RPR model includes the s-channel resonances 
Sii(1650), Pii(1710), Pi3(1720), and Pi3(1900). Apart from these, the model includes either 
a Z}i3(1900) or Pii(1900) state in the K^K channel. In detailed comparisons with the sep- 
arated structure functions sl and beam- recoil transferred polarization data from CL AS ja] , 
only the Di3(1900) assumption could be reconciled with the data, whereas the Pii(1900) 
option could clearly be rejected 25|. Note that the CLAS electroproduction data strongly 
suggest a reaction mechanism for K'^K dominated by t-channel exchange, however there are 
obvious discrepancies with the Regge predictions, indicative of s-channel contributions. 
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IV. EXPERIMENT AND DATA ANALYSIS 



A. Experimental Apparatus 

The data included in this work were taken in 1999, using the high duty factor electron 



m 



beam at Jefferson Lab and the CEBAF Large Acceptance Spectrometer (CLAS) 
Hall B. A longitudinally polarized 2.567 GeV electron beam with a current of 5 nA was 
incident upon a 5-cm-long liquid-hydrogen target with a density of 0.073 g/cm^, resulting 
in a luminosity of ~ 10^"^ cm~^s~^. The electron beam polarization was measured regularly 



throughout the experiment with a coincidence M0ller polarimeter [38|. The average beam 
polarization was measured to be 67.0±1.5%. 

CLAS is a large acceptance spectrometer used to detect multi-particle final states. Six 
superconducting coils generate a toroidal magnetic field around the target with azimuthal 
symmetry about the beam axis. The coils divide CLAS into six sectors, each functioning 
as an independent magnetic spectrometer. Each sector is instrumented with drift chambers 



(DC) to determine charged-particle trajectories 39|, scintillator counters (SC) for time-of 



flight measurements |40|, and, in the forward region, gas-filled threshold Cerenkov counters 



(CC) for electron/pion separation up to 2.5 GeV 4l| and electromagnetic calorimeters (EC) 
to identify and measure the energy of electrons and high-energy neutral particles, as well as 



to provide electron/pion separation above 2.5 GeV 42|. The trigger for the data acquisition 
readout of CLAS was a coincidence between the CC and EC in a given sector, which selected 
the electron candidates. For the data sets used in the present work, the total number of 
triggers collected was 530 M and 370 M for the two torus current settings of 1500 A and 
2250 A, respectively. These two data sets were combined together for the present analysis. 



B. Data Binning 

The data were binned in a four-dimensional space of the independent kinematic variables, 
Q^, W, cos 6'!^, and 0. Table [T] gives the binning in the variables Q^, W, and cos^l^, while 
was binned in eight, equal-sized bins running from -180° to 180°. A small fraction (< 5%) 
of the (f) bins have been excluded from this analysis due to their low acceptance in CLAS. 
A point was rejected if its acceptance was less than 2.0% (absolute) or less than 10% of the 
average acceptance over all bins at the same Q^, W, and cos 9}^. These tend to be the bins 
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adjacent to = 0° where the asymmetry is small because of the sm0 dependence seen in 
Eq.Q, and therefore their absence has little effect on the extraction of ctlt'- 
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-U.ZO 






1.75 to 1.80 


1.775 


-0 10 to 20 


0.05 






1.80 to 1.85 


1.825 


0.20 to 0.50 


0.35 






1.85 to 1.90 


1.875 


0.50 to 0.80 


0.65 






1.90 to 1.95 


1.925 


0.80 to 1.00 


0.90 






1.95 to 2.00 


1.975 










2.00 to 2.10 


2.050 







TABLE I: Ranges and centers of the kinematic bins used in this analysis. Note that the bin 
from 0.8 to 1.3 GeV^ was bin centered to the value of 1.00 GeV^ in this work (see Section [IV Hp . 



C. Particle Identification 

The p(e, e'K^)K reaction was isolated by detecting the scattered electron, e', and kaon, 
with CLAS, and reconstructing the hyperon via the missing mass technique. Electrons 
were identified by producing an electromagnetic shower in the EC accompanied by a signal in 
the CC. The electron energy deposited in the EC for all electron candidates was required to 
be consistent with the momentum measured by the track reconstruction in the DC. Electron 
and pion separation was also made by distinguishing between their different interaction 
modes in the EC. The start time of the interaction was then obtained by calculating the 
difference between the time measured by the SC and the flight time measured by the DC. 
This measured start time was combined with the hadron momentum and the path length 
measured by the DC to determine the hadron mass. 

Corrections to the electron and kaon momenta were devised to correct for reconstruction 
inaccuracies. These arise from the relative misalignments of the drift chambers in the CLAS 
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magnetic field, as well as for uncertainties in the magnetic field map employed during charged 
track reconstructions. These corrections were typically less than 1%. 

Due to the small fraction of events containing kaons in the CLAS data, a pre-selection of 
kaon events based on preliminary particle identification was made. Here the kaon candidates 
were selected by choosing positively charged particles with a reconstructed mass between 
0.3 and 0.7 GeV. Because the relative momentum resolution of CLAS becomes poorer with 
increasing momentum, a momentum- dependent mass cut was used. Fig. [3^ shows the recon- 
structed hadron mass as a function of momentum along with the cut used. Fig. [3)d shows the 
projected hadron mass distribution for all hadrons that passed the pre-selection criterion. 

3 1 

00 

o 

0.8 
0.6 
0.4 
0.2 



O 0.5 1 1.5 2 2.5 O 0.25 0.5 0.75 1 

Momentum (GeV) Mass (GeV) 

FIG. 3: (Color online) Reconstructed mass for positively charged particles. The left figure shows the 
mass plotted against the measured momentum. The lines show the mass cuts used to identify kaon 
candidates. A logarithmic yield density scale is employed. The right figure shows the reconstructed 
mass. These spectra were made from the kaon-filtered data files. The kaon peak is enhanced relative 
to the pion and proton peaks in the range from 0.3 to 0.7 GeV due to the filtering condition. 

Hyperons are identified by using the four-momenta of the electron beam, scattered elec- 
tron, and the candidate. The missing mass distribution contains a background that 
includes a continuum beneath the hyperons from multi-particle final states with misidenti- 
fied pions and protons, as well as events from ep elastic scattering (protons misidentified as 
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kaons) and events from vr+n final states (pions misidentified as kaons). The elastic events 
are kinematically correlated and show up clearly in plots of 6*j^ versus missing mass and 6k 
versus (Figs. and b, respectively). A cut on the elastic band in the 6^ (lab angle) 
versus plot removes them without a significant loss of hyperon yield. The ir^n events 
are removed with a simple missing-mass cut in which the detected hadron is assumed to be 
a pion. The resulting hyperon missing-mass distribution over the entire kinematic range is 
shown in Fig. [5l Both the A(1116) and S'^(1193) hyperons are apparent, along with several 
higher mass hyperons. 




0.6 0.8 ■ 1.2 1.4 1.6 0.5 1 1.5 



MM(eK*) (GeV) Q' (GeV") 

FIG. 4: (Color online) (a) O'^ vs. p{e,e'K^) missing mass showing ep elastic events and e'lr^n 

events. Tiie vertical bands correspond to ground state A(1116) and S''(1193) hyperons, and the 

S'^(1385)/A(1405) hyperons. (b) 6^ (lab angle) vs. for p[e,e'K^) events showing the ep elastic 

events and the cut used to remove them. 

D. Background Corrections 

To remove the multi-particle final-state background channels such as e'pvrvr, the phase 
space background was modeled by selecting the tails of the pion and proton mass distribu- 
tions. To do this, hadrons in the mass region from 0.275 to 0.725 GeV but outside of the 
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FIG. 5: (Color Online) Missing mass for p{e, e'K~^) summed over the entire range of Q^, W, cos 
and (p of the data before and after removing the elastic ep events and the e'TT~^n events. 

momentum-dependent kaon mass cuts were selected. Background missing mass distribu- 
tions were calculated for these particles assigning them the kaon mass. These background 
distributions were fit to the p{e, e'K^) missing mass distribution using a maximum log like- 
lihood method appropriate for low statistics. The fraction of each background distribution 
present in the data was thus estimated, and the normalized background contributions were 
subtracted from the data. Fig. [6] shows the missing mass distributions for two representative 
bins with the fitted background distributions overlaid. The A hyperon yields are the number 
of events in the background-subtracted missing mass spectra in the mass range from 1.095 
to 1.165 GeV. 



E. Detector Efficiency and Acceptance 

Geometric fiducial cuts were used in order to ensure that all final state charged particles 
were detected within the volume of CLAS where the detection efficiency is relatively large 
and uniform. These cuts remove the edges of the CLAS detectors and depend upon the 
momentum of the particles, as well as the torus magnetic field setting. The response of 



the CLAS detector was simulated using GSIM, a GEANT-based 43|] simulation package 
for CLAS, that combines the geometrical configuration with the inefficiencies of the various 
parts of the detector. Monte Carlo techniques were used to generate p(e, e'K^)h. events for 
each helicity state of the incident electrons by including the helicity-dependent fifth structure 
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FIG. 6: (Color Online) Examples of the background fit results for two typical kinematic bins with 
h = +1. The plots show the raw missing mass plots with the fitted background (solid histogram) 
overlaid. 



function in the Mart and Bennhold model 
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36|. 



Acceptance correction factors were obtained for each kinematic bin of Q^, cos 6*1^, 
and 0, and the two torus field settings, and the effect of the acceptance corrections on 
the helicity-dependent asymmetries was examined. In the limit of large statistics in the 
Monte Carlo simulation, the corrected asymmetries are indistinguishable from the uncor- 
rected asymmetries. One should not expect any helicity dependence to the CLAS acceptance 
outside of negligible bin-migration effects. Thus, the acceptance correction was observed to 
cancel out (within the statistical uncertainties of both the data and Monte Carlo) in the 
asymmetry, and no acceptance corrections were applied to the asymmetry measurements. 
However, a systematic uncertainty associated with not including acceptance corrections has 
been estimated (see Section IIVHI) . 

F. Radiative Corrections 

Radiative corrections were performed on the extracted reaction yields using the exact 
calculation for the exclusive approach by Afanasev et al. 4^]. This approach is based on 
the covariant procedure of infrared-divergence cancellation by Bardin and Shumeiko 45|. 
The exclusive approach is used to correct the cross section, not only in terms of the leptonic 
variables, but also the hadronic variables in exclusive electropro duct ion. These calculations 
were adapted for kaon electroproduction in this work with a cross section model that included 



18 



a contribution from ctlt'- The radiative correction factors were the ratio of the Born and the 
radiative cross sections, described in terms of four kinematic variables, Q^, W, cos 6*1^, and 
0. The radiative corrections are up to 30% for a given hehcity state but essentially cancel 
out in the asymmetries. 



G. Extraction of ctlt' 

The extraction of a^T' requires knowledge of both the asymmetry A^t' and the unpolar- 
ized cross section ctq, which can be seen by rearranging Eq.dH]) as: 

= = ctlt' sm (p. [11) 



Alt' is determined by forming the asymmetry of the K^h. yields for the positive and negative 
beam helicity states {h = ±1) as: 

da da ^ . I 

where is the cross section given by Eq.([7]) and A^^ correspond to the corrected yields for 
the positive and negative helicity states. The electron beam is partially polarized, therefore, 
the measured asymmetries are also scaled by the measured beam polarization, Pb. 

A correction for the beam charge asymmetry (differences in the integrated beam charge 
for the different helicity states of the beam) is also included. This is an extremely small 
correction and was measured to be 2.99 x 10~^. It was determined by measuring the helicity- 
dependent yield ratio for ep elastic scattering, which, outside of parity-violating effects, is 
exactly one. 

The unpolarized differential cross sections, (Tq, for the p(e, e'K^)K reaction that are used 

in this work, are the published CLAS results from the same data set l5|. The data from 
ft 

Ref. U were bin centered in Q^, W ^ and cos^l^-. In that analysis, ctq was measured with the 
same binning in the variables Q^, and cos 6*^, and the (^-dependent cross sections were 
then used to extract the structure functions, au = ctt + ^cl, arr, and ctlt- 

In order to smooth out the statistical fluctuations of the unpolarized cross section data, 
a two-dimensional simultaneous fit in and cos^^ of the data has been done. The resulting 
fitted (p- and cos ^^-dependent cross sections have been used in the extraction of aiT'- The 
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measured 0-dependent cross section in a given bin 0* is the cross section ctq averaged over 
the span of the bin from (f)\ to (ff^ (upper and lower hmits of the bin), and is given by: 

^0 = XTT / i^u + c+aLT cos + eaxT cos 20) c/0 (13) 
^0 J<pi 

^ (auA(j)' + c+aLT (sin0t - sin0j) + ^arr (sin20t - sin20j)) , 



where A0* = 0^ — 0; and c+ = a/ e(l + e). 

In addition to the trivial dependence, the unpolarized cross section has some unknown 
cos^^ dependence. It has been assumed that each of the separated structure functions can 
be described by a third-order polynomial in x = cos 9^ as: 

au = Uo + Uix + U2X^ + U3x\ (14) 
C+aLT = LTq + LTix + LT2X^ + LT3X^ (15) 
ecTTT = TTo + TTix + TTaX^ + TTgxl (16) 

Samples of the resulting fits are shown in Fig. [71 In each plot, the black solid line is the 
best fit and the dashed lines represent a ilcr error band extracted from the error matrix of 
the fit. As expected, the error band is smaller than the uncertainty of the nearby data points. 
This leads to a smaller contribution to the uncertainty of ctlt' than if the ctq data were used 
directly. The red/light dashed lines in the figures are from using the one-dimensional 
fits used in the structure function separation of Ref. [sl. The one-dimensional fits are 
very similar to the simultaneous 0/cos^^ fits and usually fall within the error band, while 
the unpolarized structure functions also agree well with those extracted in Ref. j^. This 
parameterization of the cross section is then used to determine the 0-dependent cross section 
averaged over the same bin size as each corresponding asymmetry point. 

As with the cross sections, the measured asymmetries are the average values over the 
span of the given bins. Integrating Eq.Q over the size of the bin results in: 

A Ameas slu A0 

COS (pi — COS (pu 

The asymmetry A^t' has not been corrected for the finite bin size in the variables Q^, W , 
and cos 6*1^. As will be discussed in Section TlV HI such corrections are very small compared 
to the uncertainties, and are very sensitive to the model choice. Therefore, a systematic 
uncertainty associated with not making this correction has been estimated. 
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FIG. 7: (Color online) Fits to the unpolarized cross section ctq vs. (p for our six cos O'^ points for 
VF=1.875 GeV and Q^=0.65 and 1.00 GeV^ IfoP bottom panels, respectively). The data and 
red/light short-dashed curves are from Ref. 5|]. The solid black curve in each plot is from the fit 
described in the text and the dashed curves represent a ztlo" error band around the fit. 



To extract ctlt', a simple sine fit was performed according to Eq. flTT]) . where the kinematic 
factor a/ e(l — e) has been calculated at the bin-centered value of and W for each bin (see 
Table [T]). Samples of the data and the resulting fits are shown in Fig. [81 The solid curves 
are the fit result and the dashed curves indicate the ±lo" error band from the fit. 

The error bars on the data points are a combination of the contributions from both Alt' 
and (To, and are given by: 

S{ALT'(ro) = \l{ALT'5aof + {ao5ALT'f. (18) 

The uncertainty 5 Alt' is the quadrature sum of the statistical and 0-dependent systematic 
uncertainties (see Section HVHI for details), while 5aQ comes from the fit of the cross sections 
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FIG. 8: Measured asymmetries multiplied by the unpolarized cross section and divided by the 
kinematic factor c_ = \J e(l — e) vs. 4> for our six cos0|^ points for a typical kinematic bin of 
VF=1.875 GeV and Q^=0.65 and 1.00 GeV^ (top and bottom panels, respectively). The solid 
curves show the results of the sin</> fits and the dashed lines show the ztlo" error band from the 
fits. 

described above. 

H. Systematic Uncertainties 

Various sources of systematic uncertainty that affect the measured asymmetries Aj^t' and 
the extracted structure functions oi^t' are considered in this analysis. The sources of sys- 
tematic uncertainty that affect the measured asymmetries include uncertainties due to yield 
extraction, fiducial cuts, acceptance corrections, radiative corrections, and the beam charge 
asymmetry. These are uncorrelated point-to-point uncertainties. Scale-type uncertainties 
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affect (Jlt' only and include the bin centering and beam polarization uncertainties, as well 
as the systematic uncertainties in the measurement of the unpolarized cross section, (Tq. 
Table HI] summarizes the various systematic uncertainties that affect Alt' and ct^t'- 

In the case of the 0-dependent uncertainties, all but the yield extraction uncertainty are 
dominated by statistical uncertainties. The uncertainty due to the background-subtraction 
(yield extraction) procedure was estimated to be the same as was determined in the cross 
section extraction procedure [sl- In that analysis, various changes to the procedures were 
studied, such as changing the histogram bin size in the fitting procedure and using different 
forms for the background shape (e.g. using both misidentified pions and protons, only 
misidentified pions, and only misidentified protons), and it was concluded that all systematic 
effects get larger in direct proportion to the size of the statistical uncertainty. When statistics 
are good (roughly 100 counts/bin), the residual systematic uncertainties are very small. It 
has been determined that the remaining systematic uncertainty due to the yield extraction 
is roughly equal to 25% of the size of the statistical uncertainty in any given bin (defined by 
Q^, W, and cos 6'^). This uncertainty was added linearly to the statistical uncertainty for 
the helicity-dependent yields (i.e. the overall statistical uncertainty was increased on each 
yield by a factor of 1.25). 

In order to estimate the uncertainties due to the fiducial cuts, acceptance, and radiative 
corrections, the corrected or nominal asymmetries were compared to the asymmetries that 
resulted from using either an alternative correction or cut. The RMS width of the difference 
between the nominal and alternative asymmetries, weighted by the statistical uncertainty 
of the asymmetry, was determined, and this was used as the estimate of the systematic 
uncertainty. For the acceptance effect, the difference between using no acceptance correc- 
tion (nominal) and applying an acceptance correction was studied. This is certainly an 
overestimate in this case, however, the uncertainty is small compared to the other sources 
of systematic uncertainty and much smaller than the statistical uncertainty. For the fidu- 
cial cut uncertainty, the extent of the fiducial cuts was varied over a large range. The 
resulting asymmetries were compared to the nominal asymmetries. Finally, for the radiative 
correction uncertainty, two different models were used as input to the radiative correction 
code. It has been implicitly assumed that the correction method is dominated by model 
uncertainties. 

The uncertainties for the background subtraction, fiducial cuts, acceptance, radiative 
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corrections, and beam charge asymmetry, are absolute uncertainties and were added in 
quadrature to the statistical uncertainty of each asymmetry data point before the extraction 
of a^T' ■ 



Type 


Source 


Systematic Uncertainty Alt' 


Alt' 
4> dependent 


Yield Extraction Method 


0.25xstat. uncertainty of yield 


Acceptance Function (GSIM) 


0.033 


Fiducial Cuts 


0.027 


Radiative Corrections 


0.009 


Beam Charge Asymmetry 


1.45 X 10"^ 


0"LT' 


Bin Centering 


5.8 nb/sr (absolute) 


Beam Polarization 


^cTLT' = 0.023cTLr. 


Unpolarized Cross Section 


0.124 ult', 0.115 ULT' 



TABLE II: Summary of the systematic uncertainties applied to Alt' and glt'- The two entries 
for the unpolarized cross section uncertainty are for the Q2 =o.65 and 1.00 GeV^ data sets. 



The end result of this analysis is the extraction of the fifth structure function, Glt', at 
specific points in Q^, and cos 6'^- using Eg. (ITT]) . These kinematic points are listed in 
Table [T] (note that our bin "center" for the bin from 0.8 to 1.3 GeV^ was 1.00 GeV^ as 
given by Ref. js], and not the true center at (5^=1.05 GeV^). The beam-helicity asymmetry, 
however, is sorted into particular bins of Q^, and cos 6'^, and thus a bin-centering 
correction must be considered to extract aLT' at specific kinematic points. The bin-centering 
correction would be applied to the binned asymmetries as: 

Al^> = Alt' ( ^ ) = Alt' ■ BC, (19) 

V / model 

where Af^, represents the bin-centered beam-helicity asymmetry and Alt' represents the 
bin-averaged asymmetry. To determine the bin-centering correction factors BC for this 
analysis, a model of the CLAS acceptance in Q"^ vs. W was developed to account for the 
partially filled bins. The BC factors necessarily rely on a model of Alt', where A^^^^ is 
the asymmetry calculated at a specific kinematic point (Q^, W , cos6^, 0) and A^!^, is the 
calculated bin-averaged asymmetry. The BC factors were determined using the hadrody- 
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namic model of Mart and Bennhold 



231 ] as a starting point. Within the framework of this 



model, several different choices of elementary reaction models are available with different 
ingredients, such as the resonant amplitudes included, as well as the functional forms for 
the meson and baryon form factors, i.e. the form factor, the K^K*^'-) transition form 
factor, and the A magnetic form factor. The differences between the structure functions 
derived using the different models for the bin-centering corrections on the asymmetries are 
quite small and none is clearly preferred by the asjTiimetry data. The assigned systematic 
uncertainty associated with the bin-centering corrections was chosen to be the largest RMS 
width of the (Jlt' differences using the different models. 

The relative systematic uncertainty due to the beam polarization measurement for the 
data sets used in this analysis is estimated to be ^ = 0.023. The estimated uncertainty on 
a^T' due to the systematic uncertainty on the beam polarization is given by: 



, .SPh 

>2 

b 



6aLT = \A^^r\-^ = \aM-f^- (20) 



The resulting values of aLX' also have an additional uncertainty associated with the sys- 
tematic uncertainty in the measurement of the unpolarized cross section, a^. The estimated 
systematic uncertainties for the cross sections are given in Ref. j^, which result in corre- 
sponding systematic uncertainties of 0.124(7^^/ and O.llSuir/ for Q"^ = 0.65 and 1.00 GeV^, 
respectively. The quadrature sum of the uncertainties due to the bin-centering correction, 
beam polarization, and unpolarized cross section, are shown by the shaded bars on the 
results (see Section IV] Figs. [9]to [T2|) . 



V. RESULTS AND DISCUSSIONS 

The angular dependence of aLT' for various W points for the two points is shown 
in Figs. [9] and [TOl along with comparisons to several model calculations. The lower 
data shown in Fig. [9] is rather fiat over the full range of energy and angle, with no strong 
structures visible. Unlike the low data, a strong W and angular dependence is observed 
in the higher data (Fig. [TOl) . The angular dependence shows an interesting peaking at 
middle angles for the lowest W point {W=1.Q5 GeV), while a rapid sign change is seen at 
both W= 1.875 and 1.925 GeV at central angles. 

The extracted aLT' structure function results are shown as a function of W for various 
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cos^^ points for the points at 0.65 and 1.00 GeV^ in the top panels of Figs. [TT]and[T2l 
respectively. For the lower data, Fig. [11] shows that for the four backward-most kaon 
center-of-mass scattering angles (cos =-0.60, -0.25, 0.05, and 0.35), (7^^/ exhibits a smooth 
energy dependence with a fall off of the structure function to zero at the highest W points. 
In the forward kaon scattering angles, cos6'^=0.65 and 0.90, where the reaction is expected 
to be dominated by t-channel exchange, ctlt' is consistent with zero to within the rather 
large error bars of the data, and no obvious structures are present. This might indicate the 
dominance of a single t-channel exchange. 

For the higher data (see Fig. [T2l) . the range of W is limited by the CLAS acceptance. 
Here the W dependence of Oi^j^i is similar to the lower data at the more forward angles, 
cos6'^=0.65 and 0.90. However, there is a notable feature in the W dependence in the back- 
ward and middle kaon angles. At cos ^|^=-0.25, 0.05 and 0.35, the data show an interesting 
interference feature around 1.9 GeV, with a rapid change of sign at cos6'^=0.05 and 0.35. 
While at the very backward angles, cos6'|j-=-0.60, a strong enhancement is seen at about 
W=1.7 GeV with a flat response for higher W. For both of the values, aiT' goes to zero 
at higher W. 

The results are compared with calculations from the MB isobar model 46| (blue long 
dashes), the GLV Regge model 47| (red dash-dot), the Ghent RPR model |48| including a 



-Di3(1900) state (green sohd), and the Ghent RPR model |48|] including a Pii(1900) state 
(green short dashes). In general, none of the available models fully describes these data 
over the Q^, W, and cos^^ ranges measured. The MB and GLV models under predict 
the strength of aiT', although they qualitatively follow the trends of the data. From the 
comparisons, the RPR model including the Pii(1900) state is clearly ruled out as already 
indicated in Ref. {25!, but the RPR model including the -Di3(1900) state seems to best 
describe the data qualitatively. However, each of these models misses key features of the 
data. The disagreements with the isobar models (MB and RPR) may not be too surprising 
as they have not been fit to these data. Therefore, the aiT' structure function provides for 
additional new constraints on the model parameters. 

A direct comparison of the measured polarized structure function aiT' with a^x from 
Ref. [5] can reveal some interesting features of the data. This is shown in Figs. [TT] and [T2l 
where the polarized structure function aLT' is plotted as a function of W for various cos 6"^ 
bins and compared with at the same kinematic points. The magnitudes of the two 
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FIG. 9: (Color online) Polarized structure function a^T' (nb/sr) vs. cos^^ for 0^=0.65 GeV^ 
and W points as indicated. The curves are calculations from the MB isobar model 46|] (blue long 
dashes^the GLV Regge model \^\ (red dash-dot), the Ghent RPR model inclu ding a Di3(1900) 



state 



(green solid), and the Ghent RPR model including a Pii(1900) state 



48( 1 (green short 



dashes). The shaded bars indicate the estimated overall systematic uncertainty on the results. 



structure functions are comparable in both the lower and higher Q"^ data, although a^T' has 
larger uncertainties. In the lower data at the most backward kaon center-of-mass angle, 
cos ^^=-0.60, (Tlt is essentially zero, while ctlt' is clearly non-zero. At cos ^^=-0. 25, 0.05, 
and 0.35, aLx is similar in shape and magnitude to o"lt', but with an opposite sign. At the 
very forward kaon center-of-mass angles, cos6'^=0.65 and 0.90, aiT' is consistent with zero, 
while a^T is non-zero. In the higher data, and for backward and middle kaon scattering 
angles, aiT' has some significant deviations from a smooth behavior, indicating significant 
interferences. However, the shapes of ctlt' and a^x are quite different. 

All of the data included in this work have been entered into the CLAS physics 
database [49|]. 
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FIG. 10: (Color online) Polarized structure function (Jlt' (nb/sr) vs. cos^l^- for Q^=im GeV^ 
and W bins as indicated. The curves are as indicated in Fig. [9l 

VI. CONCLUSIONS 

The first measurements of the structure function air' for K^A electroproduction have 
been reported. The data span a range in W from threshold to 2.05 GeV for two points at 
0.65 and 1.00 GeV^, and span nearly the full center-of-mass angular range of the final state 
In this analysis, the energy and angular dependence of ctlt' have been investigated. 
aiT' is found to be comparable in size to the unpolarized cross sections. The structure 
function is surprisingly featureless with energy and angle for the lower data, while the 
higher data indicate rather strong interference affects near threshold and at ly of 1.9 GeV 
for central angles. a^T' is consistent with zero at more forward angles and at higher values 

of ly. 

The data have been compared with several different model calculations. The GLV Regge 
calculation generally under predicts the data. This is perhaps not too surprising given that it 
includes no explicit s-channel processes, which are expected to show clear signatures in a^T'- 
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FIG. 11: (Color online) Polarized structure function (y^^l as a function of W for various cosO*^ 
(upper panel) compared to the measured aiT from Ref. 5| (lower panel) for (5^=0.65 GeV^. Plots 
are shown for cos 0|^=-O.6O, -0.25, 0.05, 0.35, 0.65, and 0.90 (left to right). The curves are as 
indicated in Fig. [21 

Comparisons to the MB isobar model and RPR hybrid isobar/Regge model (which did not 
include any CLAS electroproduction data in their fits) indicate that the model parameters 
need to be tuned in order to reproduce the overall average strength seen in a^T'- A bigger 
challenge for these models is to explain the strong interference signatures in the data. Even 
though the CLAS aiT and aiT' data have rather sizable statistical uncertainties, the data 
do have a good deal of discriminating power with regard to certain assumptions about which 
resonant states are included. 

The ctlt' results were also compared with the results of the measurements for 
at the same kinematic points. While ctlt' has larger uncertainties, the magnitudes of the 
two structure functions are comparable. In our lower data, a^T' is quite smooth with 
W and cos^lf. However, at the high value, ctlt' shows strong interference and/or FSI 
signatures at middle and backward kaon scattering angles. Together, these two observables 
provide more complete information on the amplitudes underlying the longitudinal-transverse 
response for this reaction. 

The question of the presence of any new resonances must wait for further work with the 
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FIG. 12: (Color online) Polarized structure function oirri as a function of W for various cos9'^ 
(upper panel) compared to the measured aiT from Ref. 5| (lower panel) for Q^=1.00 GeV^. Plots 
are shown for cos 0|^=-O.6O, -0.25, 0.05, 0.35, 0.65, and 0.90 (left to right). The curves are as 
indicated in Fig. [21 



existing hadro dynamic models and partial wave analyses applied to the full range of the 
data. Fortunately, the new information presented here will impose reasonable constraints 
on the amplitudes used to describe electroproduction of -R''^A and K'^T,'^ final states, making 
these models more reliable for future interpretation and prediction. 
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